Quark matter under strong magnetic field in chiral models. 
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The chiral model is used to describe quark matter under strong magnetic fields and compared 
to other models, the MIT bag model and the two flavor Nambu-Jona-Lasinio model. The effect of 
vacuum corrections due to the magnetic field is discussed. It is shown that if the magnetic field 
vacuum corrections are not taken into account explicitly the parameters of the models should be 
fitted to low density meson properties in the presence of the magnetic field. 
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I. INTRODUCTION 



Magnetars, neutron stars with very strong magnetic fields of the order of 10 14 — 10 15 G at the surface, are sources 
of very energetic electromagnetic radiation, mainly gamma and X rays [1-3]. Presently, about 21 of these objects 
have been detected, most of them as soft gamma repeaters (SGRs) and anomalous X-ray pulsars (AXPs) [4]. 

It has been argued [5] that strange quark matter (SQM), i.e. quark matter with strangeness per baryon of the order 
of unity, may be the true ground state of hadronic matter. This could imply that compact stars are mainly quark 
stars (see also [6]). Magnetars as compact quark stars have been first investigated in [7], where the MIT bag model 
[8] was applied to obtain the equation of state (EOS) of stellar quark matter under magnetic fields as strong as 10 18 
G. The EOS for magnetized quark stars described within the MIT bag model and taking into account the anomalous 
magnetic moment of quarks (AMM) was studied in [9]. The su(2) version of the Nambu-Jona-Lasinio model (NJL) 
, [10], an effective model which includes the chiral symmetry, was applied to study the stability of quark matter under 
CNJ ■ very strong magnetic fields in [11]. The same model and its su(3) extension were used to describe quark stars with 
very large magnetic fields in [12, 13]. In these papers, the phenomenon of magnetic catalysis within the NJL models 
has been discussed. Magnetic catalysis is one of the most important effects of the magnetic field in quark models 
with chiral symmetry and corresponds to the enhancement of the chiral symmetry breaking in the magnetic field [14] . 
Another nontrivial effect of the magnetic field is the possibility that strong magnetic fields can turn a crossover into 
a first order QCD transition [15]. 

A different quark model with chiral symmetry, the chiral model of pions and quarks or nucleons [16-18], also known 
as linear sigma model, was applied to study the high density npe/x matter with 7r°-condensation [19]. In the present 
work we will use the same model to describe both symmetric quark matter and stellar quark matter under strong 
magnetic fields. No pion condensation will be considered. The results will be compared with the MIT and NJL 
models. In particular we will investigate the inclusion of vacuum corrections due to the strong magnetic field. 

The linear sigma model coupled to quarks and the Polyakov loop has been recently used to study the phase diagram 
of hot QCD in a strong external magnetic field [20]. There it was shown that a strong magnetic field could give rise 
to a splitting of the dcconfincment and chiral transitions if the B dependent vacuum corrections were included. These 
results agree well with the diagram coming from the NJL model [21], where the vacuum corrections are present 
authomatically, and with the results of the lattice calculations [22]. 

The MIT bag model describes quarks as a free gas of quarks already in a chiral restored state. The bag pressure 
provides the confinement and is just a parameter which can be fixed from the nucleon sector. The problem of chiral 
symmetry restoration is beyond the scope of this model. Both NJL and the chiral model are described by chiral 
symmetric Lagrangian densities and a vacuum state with spontaneously broken chiral symmetry. In the chiral model, 
the chiral condensate plays the role of the bag pressure and its value in the vacuum is fixed from the pion decay 
constant, which is well known, and the sigma mass. The connection between the MIT bag model and the chiral model 
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has been discussed in [23]. In the NJL model the model parameters are fixed by fitting the pion decay constant and 
the quark condensate. 

In section II we make a brief review of the three models and their corresponding EOS under the effect of a magnetic 
field and discuss how the model parameters are fixed for a finite magnetic field. Results are discussed in section III 
and conclusions are drawn in section IV. 



II. QUARK MODELS 

In the present section we give a brief review of the quark models including the effect of a strong magnetic field used 
in this study: the chiral sigma model, the su(2) NJL model and the MIT bag model. 



A. Chiral sigma model 

We consider the chiral sigma model for quarks interacting with an external magnetic field. The chiral symmetric 
Lagrangian density reads [16, 18] 

C = i> } {i-f^ - qipAi* - g (a + i~fyT.it)) ipf 

+ \d^a + ^n°d^° + £> M 7r + £P7r- - U (a, n) - \f^F„ v , (1) 

where ipf is the quark field, a and 7r°, tt = (it 1 ± i-K 2 )/^/2 are the meson fields, = + ieA^, g is the quark-meson 
coupling constant, q = (3 + T 3 )e/2 and = (0, 0, Bx, 0) refers to an external magnetic field along the z-axis. 

The potential functional U is a " Mexican hat" potential, which leads to spontaneous chiral symmetry breaking and 
is included to reproduce the vacuum expectation value of the sigma field. For exact chiral symmetry (m n = 0), the 
potential is 

U(a,n) = ^(a 2 +n*-f n ) 2 . (2) 

where A 2 = — ^, and m a is the mass of the cr-meson. The vacuum expectation value of the field a is (a) = f v , 

where f v = 93 MeV is the pion decay constant in the absence of a magnetic field. In the normal phase with no pion 
condensation, the energy density is given by 

£ = E -^TT- E a " / dpJp}+m> + 2v\q f \B + U(cT,0) (3) 

f=u,d v=V 70 

1 S 

where m = ga, and v = n+ — — sgn(q) - — 0, 1, 2, . . . enumerates the Landau levels (LL) of the fermions with electric 
charge q, the factor a v — 1, (2) for v = (y > 1) takes care of singly degenerate zeroth LL and doubly degenerate 



LL levels with u > 1 and the Fermi momentum of LL v is PfJv = \J V/ ~ Mf(v, B). The coefficient N c — 3 stands 

for the color degeneracy. We use m a = 1200MeV. 

The energy minimization for each baryon density with respect to the a field, gives the so called gap equation 

v NMAB ^fr d Pz +!E = n m 

kU 2 " 2 » Vpl+™* + 2v\q f \B + da U 

The mass in the vacuum is m = gf n . In the presence of a strong magnetic field the vacuum properties are affected 
and a quark vacuum correction should be included [12, 20]. This correction was not considered in [19]. Since the 
magnetic fields discussed were below 6 x 10 18 G, we do not expect a large effect coming from this term. We will take 
into account this correction in two different ways: a) we will redefine the constant and suppose all vacuum effects 
are described by the "Mexican hat" , including the magnetic field contribution; b) we consider that the "Mexican hat" 
potential does not include the magnetic field vacuum contribution and we will add this contribution explicitly as an 
extra term to (2) just like it was done in [20]. 
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B. The MIT bag model 

In the presence of a strong magnetic field, the energy density and quark density within the MIT bag model, are 
given by [24] 



X- N c \q f \B x ^ 

2^ 4^2 2^ a " 

f=u,d v=0 



HfPFjv + M 2 f (v, B) In 



M f (u,B) 



+ Bag, 



/=«>' 



X ^ N c \q f \B ^ 
Pi = 1^ 2 7T 2 ^ a ^ F ' /!y 

y=0 



(5) 
(6) 



where Mf(v, B) = + 2v\qf\B and v runs over the allowed LL, m q is the quark mass and Bag represents the bag 

pressure. We only consider flavors u and d. 



C. The su(2) NJL model 



We consider the two flavor NJL model defined by the following Lagrangian density [12, 13] 

c = c f - \f, v f^ 

where the quark sector is described by the Nambu-Jona-Lasinio model 

+ G [(^fPf) 2 + (TpfijUTTpf) 2 ] , 



(7) 



(8) 



m c = m u ~ m,d are the quark current masses; and = d^A v — d v A^ are used to account for the external 
magnetic field. Since we are interested in a static and constant magnetic field in the z direction, A^ = S^xiB. 
The energy density is given by 



£(v f ,B) = -P N + J2HPf , 
f 

where /if is the chemical potential of flavor / and the pressure is P N = P(Hf)\M(n f ) — P (0)|m(o) with 

P = On + e d - G(<j> u + <i> d f , 

For a given flavor, the Of term is given by 

d 4 p 



"1 = -?'] (2»)« 
and the condensate 4>f = (•tpfipf), so that [12, 13] 

p ^pvac _|_ pmag _|_ pmed^ 

where the vacuum contribution reads 



N c N f 
8tt 2 



M In 



A+_fA 

M 



6 A A(A 2 + e l) 



with e\ = VA 2 + M 2 , A representing a non covariant ultra violet cut off; the finite magnetic contribution is 



pmag _ \^ 



f—u,d 



Nc{\q f \Bf 
2tt 2 



1 

Q'{-hx f )- ~{x) - Xf )]nx f + -£ 



(9) 



(10) 



(11) 



(12) 



(13) 



(14) 
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with Xf — M 2 /(2\qf\B) and £'(— l,Xf) = d((z,Xf)/dz\ z =-\ where ((z,Xf) is the Riemann-Hurwitz zeta function, 
and the medium contribution can be written as 



P med sr \if\ BNc sr n 
r M - 2^ 4^-2 2^ a ^ 

f=u,d v=0 



VfPFjv - B) 2 ln 



A 4 / "I" PfJv 



M f {y,B) 



-2G(<t> u + <j> d ) 2 , 



(15) 



where Mf{v, B) = \J M 2 + 2\qj\Bv and Pf,/i/ = yV/ — ^KfO^i ^) 2 - The upper Landau level (or the nearest integer) 
is defined by 



M 2 -M 2 » 2 



2|«/|S 2| g/ |S- 
The effective quark masses can be obtained self consistently from 

M = m c - 2G{4> U + fa), 

where the condensates 4>f are given by 



A I 



with 



<j) v f ac = 



MNc 
' 2ir 2 



(16) 



(17) 



(18) 



(19) 



M\q f \BN c 

2tt 2 



]nT(x f ) - \ ln(27r) +Xf~\ (2x f - l)Hx f ) 



and 



M\q f \BN c 
2^ 2 



E 

L/=0 



ttj, In 



/ Hf+PF,fu 

\ M f (v,B) 



(20) 



(21) 



The density, p/, corresponding to each different flavor, is given by Eq. (6). 

For the su(2) NJL model we use the same parametrization as given in [25]: A = 587.9 MeV, m c — 5.6 MeV, 
m e = 0.511 MeV, — 105.66 MeV, and GA 2 = 2.44 which gives a quark vacuum mass equal to 400 MeV in the 
absence of a magnetic field. 



D. Fixing the model parameters at finite B 

In order to compare the MIT bag model with the other two models we chose Bag = 0.9 far 4 = 177.3 Me V/fm 3 
which is in between the values we get for the chiral model and the NJL model in the chiral symmetric density region. 
We note that in MFT, when the chiral symmetry is restored, i.e. ((a) — and (n) = 0), the potential functional 

A 2 

reduces to a constant energy density equal to —f^, an d the constituent quark mass vanishes, leaving free massless 

quarks. The chiral Lagrangian, in the region where the chiral symmetry is restored, can be identified by a MIT bag 

A 2 

quark model with a bag pressure Bag= ~ft- 

In order to be able to establish a comparison between the chiral model and the NJL model, we include, in the chiral 
model, vacuum corrections due to the magnetic field in two different ways: (a) by renormalizing the pion constant 
f' v constant so that the vacuum quark mass mo = g in the chiral model coincides with the one obtained with the 
NJL model. We will call this model the chiral model I; (b) by including in the potential (2) the explicit field vacuum 
contribution, U mag (a,B) = -P ma 9 7 w here P mag was defined in (14) with M = ga as was done in [20]. This will be 
the chiral model II. In this case we have for the potential 

U (a, *,B) = ^ (a 2 + tt 2 - ,/ 2 ) 2 + U mag (a, B) - U(a , 0, B) (22) 
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B=0 




FIG. 1: (Color online) The quark mass as a function of baryon number density for isospin-symmetric quark matter, and for 
several values of the magnetic field. The lines correspond to the su(2) NJL (black lines), the chiral model I with a renormalized 
f-n (blue lines) and the chiral model II with B dependent vacuum corrections defined in (23) (red lines). 



where 



U mag (a,B) = 



f=u,d 



2tt 2 



C'(-M/) - \(*t ~ Xf)lnxf + ^ 



(23) 



with Xf = (ga) 2 / (2\qf\B) as in (14), and CTo is the a field in the vacuum. The last term in (22) insures that the 
pressure goes to zero at zero density. Replacing this expression for the potential in the gap equation we get: 



E 

f=u,d 



2tt 2 



E o 

\q f \BN c 
2tt 2 



dp z 



y/pl +m 2 + 2v\q f \B gm 
^T(x f ) - i ln(27r) +*f-\ (2xf - 1) \n(x f ) 







(24) 



In Fig. 1 we compare the quark masses in isospin-symmetric quark matter as a function of the baryon density 
obtained with both models and using the two approaches described above to take into account the B dependent 
vacuum corrections. By construction the curves for the NJL and the chiral model I start at the same mass at zero 
density and decrease with density reaching the value of the current mass of the quarks in the chiral symmetric phase. 
For the chiral model II the quark vaccum mass also increases with B but not so strongly. As discussed in [11-13], in 
the NJL model, the magnetic field shifts the chiral symmetry restoration to larger densities. The same occurs with the 
chiral model II. We also note that within this model the chiral symmetry restoration occurs slower than in the NJL 
model, a feature that is already present for zero magnetic field. In the chiral model I magnetic catalysis is not so clear. 
There are two competing effects which can be identified from the gap equation: for a larger i.e. larger vacuum 
quark mass, the restoration of chiral symmetry occurs at smaller densities if no Landau quantization is present in the 
quark quasi-particle (QP) energy. This is the effect that is observed for the smaller magnetic fields. However, Landau 
quantization reduces the QP energy, and, therefore, for very strong fields, the chiral symmetry restoration occurs at 
larger densities when B increases. 
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FIG. 2: (Color online) Quark matter energy per baryon number as a function of the baryon number density for isospin- 
symmetric quark matter, and for several values of the magnetic field. The different lines correspond to the NJL (dashed lines), 
chiral model I (full lines), chiral model II (dash-dotted lines) and MIT bag model (dotted lines). 



III. RESULTS 



The main objective of the present section is to compare the different quark models and to discuss how the vacuum 
corrections due to the presence of a strong magnetic field may be taken into account in the MIT bag model and in 
the chiral model. We will first compare the properties of symmetric baryonic matter and in a second subsection, we 
will discuss the implications in the EOS of stellar quark matter obtained with both chiral models considered. 



A. Symmetric quark matter 

In Fig. 2 the energy per baryon calculated with the chiral model, NJL model, and MIT bag model, is shown for 
several values of the magnetic field. In the absence of an external magnetic field the three models give similar results 
at high densities, when a chiral symmetry restored state is the ground state of the system. Since, both the chiral 
model and the NJL model use low density meson properties to fix the parameters of the models they behave in a 
similar way. Below B = 10 19 G the magnetic field has no noticeable effect on the energy per particle. However, above 
B = 10 19 G the effect becomes stronger and the MIT gives the lowest energy per particle. This is due to the fact that 
the bag pressure, which describes the vacuum effects, was kept constant, independent of the magnetic field, and, due 
to the Landau quantization, the contribution of the kinetic energy is strongly reduced. The minimum of the energy 
per particle is shifted to higher densities but not so strongly as in the other models. In the chiral model I we see that: 
a) in the chiral symmetry broken phase the chiral model has a behavior similar to the NJL for large values of B, and 
at zero density the energy per particle of both models always coincide, b) although B dependent vaccum corrections 
are partially taken into account and the energy per particle does not decrease so much as in the MIT bag model, in 
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FIG. 3: (Color online) The bag pressure as a function of baryon number density for isospin-symmetric quark matter. The 
different lines correspond to the NJL (dashed lines), chiral model I (full lines), chiral model II (dash-dotted lines) and MIT bag 
model (dotted lines). Please notice that the y-axis scale of the bottom figures is different from the one of the top figures. 



the chiral restored phase the energy per particle is still much smaller than the predictions of both the NJL and the 
chiral model with the vacuum correction (23) for B > 5 x 10 19 G. In chiral model II, where the B depedent vacuum 
corrections are properly taken into account in the chiral model we conclude that: a) in the chiral symmetry broken 
phase the chiral model predicts for B < 3 x 10 19 G a larger energy per particle than NJL. However, for larger values 
of B the opposite occurs because the vacuum mass in this model does not increase so much with B than in the NJL 
model; b) in the chiral restored phase the energy per particle is larger in the chiral model until a very large magnetic 
field (10 20 G). 

Let us define, in all models, an effective bag pressure which corresponds to the total energy minus the kinetic 
energy contribution. In Fig. 3 we plot the effective bag pressure for the three quark models, and several magnetic 
field intensities. We have considered both approachs to the B dependent vacuum corrections in the chiral model. For 
B = 0, all models are supposed to describe the same physics at high densities which corresponds to the chiral restored 
phase: the NJL has the smallest bag pressure while the chiral model the largest one. For the MIT, as discussed 
before, we have chosen an intermediate value. For finite values of B, we conclude that: a) in the chiral symmetry 
broken phase the chiral model II has the smallest effective bag pressure although close to the one of the NJL model. 
The chiral model I gets closer to the NJL model as B increases, however, the chiral symmetry restoration occurs at 
too low densities compared with the other two chiral models; b) in the chiral symmetry restored phase the effective 
bag increases with B in all models except in the MIT bag model. Within the NJL, the effective bag increases faster 
with B than in the chiral model II and at B = 10 20 G both models coincide. In the chiral model I the effective bag 
increases too slowly and above 4 x 10 19 G it is already smaller than the one obtained in the NJL model. However, 
we should point out that the large differences occur only for very intense magnetic fields, above 5 x 10 19 G. 

The behavior of the bag pressure in the different models reflects itself on the pressure of the system. In Fig. 4 we 
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FIG. 4: (Color online) Pressure as a function of baryon number density for isospin-symmetric quark matter, and for several 
values of the magnetic field. The different lines correspond to the NJL (dashed lines), chiral model I (full lines), chiral model 
II (dash-dotted lines) and MIT bag model (dotted lines). 

plot the pressure of the gas of symmetric quark matter as a function of density for different magnetic field intensities. 
Below 10 19 G all models give results in the chiral restored phase similar to the ones obtained at B = 0: the model 
with the smallest effective bag pressure (the NJL model) has the largest pressure. Above 3 x 10 19 G the different 
effective bag pressures of the three models give rise to quite different pressures at high densities with the MIT bag 
model predicting the largest one and the chiral model II the smallest one. 

B. Stellar quark matter 

In a quark star, we must impose both /3-equilibrium and charge neutrality. The relations between the chemical 
potentials of different particles are given by 



Hd = Pu+Pe, /Xe = /V 

In terms of the neutron and the electron chemical potentials fi n and /x e , one has 



1 2 



1 1 



Pu — 3^ e ' ^ d 3^" 3^ e 



(25) 



(26) 



For the charge neutrality we impose 



Pe + P M = 3 ( 2 P« - Pd) ■ 



(27) 
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FIG. 5: (Color online) Pressure (top) and particle fractions (bottom) as a function of baryon number density for asymmetric 
quark matter, and for B — 0, 5 x 10 18 , 10 19 G. The thin line is for the su(2) NJL model and the thick line is for the chiral 
model. The zero-axis is shown in the pressure plots with a thin dotted line. In the bottom plots the black/thick lines are for 
the su(2) NJL model and the red/thin lines for the chiral model I. 



In Fig. 5, we plot the pressure and the particle fraction as a function of the baryon density obtained with the 
chiral model and the NJL model for different values of the magnetic field intensity. For the chiral model, we will 
only consider the B dependent vacuum effects through the rcnormalization of the pion decay constant, chiral model 
I. Since we will not go beyond B = 10 19 G both prescriptions introduced in section II D give similar results. 

At the surface of the quark star, defined by a zero pressure, the density is finite. The chiral model predicts larger 
baryon densities at the surface and a softer EOS, for magnetic field below 10 19 G. For larger fields the opposite 
occurs, e.g. NJL predicts larger densities at the surface. However, according to the scalar virial theorem [26] the 
interior magnetic field strength could be as large as B ~ 1 — 3 x 10 18 G so, in principle, fields stronger than the ones 
represented in Fig. 5 will not occur in the interior of compact stars. 

Considering the particle fractions obtained within both models, it is seen that the NJL model predicts a larger 
w-quark fraction, and consequently, larger electron and muon fractions. At B — 10 19 G the main effects due to the 
magnetic field occur below p = 0.25 fm~ 3 . The irregularity of the curves is due to the filling of the Landau levels, 
which may give rise to strong fluctuations on the particle fractions. There is a large increase of the it-quark, and 
correspondingly of the electron fraction, with the increase of B. It is even observed that the u quark fraction is larger 
than the d quark fraction in a small range of densities. However, we should point out that according to the pressure 
plots the density at the surface of the star will be above 0.5 fm -3 . Therefore, the strong effects on the particle fractions 
below that density will not affect the star properties, except if the star has a crust as discussed in [27, 28]. In this 
case strong effects could occur in the star crust. The onset of the muon is sensitive to the field and reflects the filling 
of the lepton Landau levels; larger muon fractions are attained at high density due to the larger u quark fractions. 
For the chiral model the muon fraction is always below 0.001 in the absence of a magnetic field. This changes for a 
magnetic field stronger than 10 18 G for densities above 1 fm -3 as seen in Fig. 5 and 6. 

In order to better understand the effect of the magnetic field on the particle fractions we plot in Fig. 6 the particle 
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FIG. 6: (Color online) Quark and lepton fractions as a function of the magnetic field for asymmetric quark matter, and for 
several values of the baryonic density (p = 0.001, 0.01, 0.3, 0.5, 0.8, 1.2 fm -3 ). The black/thick lines are for the su(2) NJL 
model and red/thin lines for the chiral model. 



fractions as a function of the magnetic field intensity for six representative densities: the densities below 0.3 fm -3 
would only occur in the crust of the star, in the case it exists, ~ 0.5 fm -3 is the surface density, and 0.8 and 1.2 
£m~ 3 are baryon densities in the interior of the star. At the surface the electron and muon fractions are larger for 
the NJL model. This has implications in the possible existence of a crust [27]. A larger electron fraction will be able 
to support a larger crust. In [29] it was shown that NJL would predict a larger electron fraction at the surface than 
the MIT bag model. The chiral model seems to behave more like the MIT bag model. At the surface of a star with 
no crust, the effect of the magnetic field starts to be non-negligible for B > 7 x 10 17 G with a clear increase of the 
electron fraction in the chiral model. This effect will occur for much smaller fields in a star with a crust as can be 
seen in Fig. 6a). As a result, it is expected that the structure of the crust will be deeply influenced by the presence 
of a very strong magnetic field. 

The effect of the magnetic field on the muon fraction within the chiral model is also clearly seen in this figure: for 
p = 1.2 fm~ 3 the muon fractions rises above 0.001 for fields a bit below 10 18 G. 



IV. CONCLUSIONS 



In the present work we have compared the properties of quark matter under a strong magnetic field described using 
three different models: the MIT bag model [8] which describes quark matter in a chiral restored phase, the chiral 
model [17, 18], and the two flavor NJL model [10] both described by a chiral symmetric Lagrangian density. We have 
discussed the effect of the magnetic field on the vacuum properties and how the parameters of the models related to 
the vacuum should be chosen in order to take into account vacuum corrections. However, these corrections are only 
important for very strong magnetic fields B > 3 x 10 19 G, which are not expected to be found in compact stars but 
could be formed as short-lived magnetic fields in relativistic heavy ion collisions playing an important role in possible 
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experimental signatures of strong CP violation and the phenomenon of the chiral magnetic effect [30]. Estimations 
done in [31] show that for the LHC energy it could be possible to get eB <~ 15 m 2 which corresponds to a field 
B~5x 10 19 G. 

It was shown that if the schematic MIT quark model is used to describe quark matter under strong magnetic fields 
the value of the bag pressure should be adjusted in order to account for the magnetic field vacuum corrections. In 
the chiral model the vacuum corrections may be taken into account by fitting the parameters of the model to a quark 
vacuum mass that includes these corrections, chiral model I, or by including B dependent vacuum corrections as 
discussed in [20], chiral model II. It was shown that for very strong fields, above 5 x 10 19 G, the chiral model I failed 
to include adequatly the B dependent vaccum correction. However, below those extreme magnetic fields chiral model 
I gave reasonable results. Although chiral model II behaves like the NJL model there are some important differences: 
its vacuum mass does not increase so fast as the one of the NJL when B increases. As a result, in the chiral model II 
the energy per particle in the chiral symmetric broken phase becomes smaller than the one of NJL for very large B 
values. In the chiral symmetric phase the effective bag pressure increases faster with B in the NJL model. Therefore, 
although the bag pressure is smaller in the NJL for B = 0, at B = 10 20 G, NJL and the chiral model II have the 
same bag pressure. 

Finally, we have applied the chiral model and the two flavor NJL model to the description of stellar matter. It was 
shown that within the chiral model quark stars will have a larger baryon density at the surface. This feature could be 
reversed for a strong magnetic field larger than 3 x 10 19 . However, so strong magnetic fields are not expected to exist 
in compact stars. Another difference is the larger (smaller) u (d) quark fractions in the NJL model. As a consequence 
the chiral model predicts much smaller lepton fractions. One of the main effects of the magnetic field is to increase 
the u quark fractions at low densities due to their larger absolute charges, and, therefore, also the electron fractions. 
This could have strong effects on the strucutre of the crust of a quark star like the one predicted in [28] . 
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